Elastic constants and ultrasonic attenuation in the cone state of the frustrated 

antiferromagnet CS2CUCI4 



Andreas Kreisel and Peter Kopietz 
Institut fiir Theoretische Physik, Universitdt Frankfurt, 
Max-von-Laue Strasse 1, 60438 Frankfurt, Germany 

Pham Thanh Cong, Bernd Wolf, and Michael Lang 
Physikalisches Institut, Universitdt Frankfurt, Max-von-Laue Strasse 1, 60438 Frankfurt, Germany 

(Dated: July 15, 2011) 

In an external magnetic field perpendicular to the plane of the layers, the quasi two-dimensional 
frustrated antiferromagnet CS2CUCI4 exhibits a magnetically ordered "cone state" at low temper- 
atures. In this state the component of the magnetic moments in field direction is finite, while 
their projection onto the plane of the layers forms a spiral. We present both theoretical and ex- 
perimental results for the magnetic-field dependence of the elastic constants and the ultrasonic 
attenuation rate in the cone state. Our theoretical analysis is based on the usual spin-wave expan- 
sion around the classical ground state of a Heisenberg model on an anisotropic triangular lattice 
with Dzyaloshinskii-Moriya interactions. Magnon-phonon interactions are modeled by expanding 
the exchange interactions up to second order in powers of the phonon coordinates. As long as the 
external magnetic field is not too close to the critical field where the cone state becomes unstable, 
we obtain reasonable agreement between theory and experiment, suggesting that at least in this 
regime magnons behave as well-defined quasiparticles. We also show that the assumption of well- 
defined magnons implies that at long wavelengths the ultrasonic attenuation rate in the cone state 
of CS2CUCI4 is proportional to the fourth power of the phonon momentum. 

PACS numbers: 43.35. -fd, 75.10.Jm, 75.30.Ds, 72.55.-|-s 



I. INTRODUCTION 

In recent years a large amount of theoretical and ex- 
perimental work has been devoted to low-dimensional 
quantum magnets on frustrated lattices!^ In these sys- 
tems enhanced quantum fluctuations and competing in- 
teractions can stabilize magnetically disordered ground 
states in certain parameter regimes. The phase transi- 
tions between ordered and disordered phases as a func- 
tion of some nonthermal control parameter, such as an 
external magnetic field, have attracted a lot of attention. 
An important realization of such a frustrated magnet ex- 
hibiting a rather complex phase diagram, including mag- 
netically ordered and disordered phases, is the magnetic 
insulator CS2CUCI4. This materi al w as synthesized and 
characterized about 15 years ago.l^ Coldea et alJ^ de- 
termined the magnetic structure and pointed out that 
it is a realization of a triangular lattice antiferromagnet 
with two different exchange couplings J and J' whose 
magnitude is less than 1 meV. The effective (7-factor 
was obtained by electron spin resonance measurement^. 
Later on, an effective model Hamiltonian was deduced 
from high field measurements.^ The observation of Bose- 
Einstein condensation of magnonJ3 and investigations on 
the phase diagram followed. An ordered cone state and 
a ferromagnetic phase which are connected by a quan- 
tum critical point were discovered. The challenge to gain 
a quantitative microscopic understanding of the exper- 
imental observations has motivated a large number of 
recent theoretical works.'^'^ 

There are two complementary theoretical points of 



view: On the one hand, it has been argued that the differ- 
ence in the two competing exchange couplings in the lay- 
ers is sufficiently large to view the system as a collection 
of weakly coupled spin chains.'^'*'*"'^''' The properties of the 
decoupled one-dimensional spin chains can then be ob- 
tained nonperturbatively using bosonization techniques, 
while the weak coupling between the chains is usually in- 
cluded using some kind of mean-field approximation. Be- 
cause in this approach one assumes that the elementary 
excitations of the system can be connected to the spinous 
of the one-dimensional spin chains, this approach is most 
natural to describe the spin-liquid sector of the phase di- 
agram. On the other hand, for the description of the part 
of the phase diagram which exhibits magnetic order, it 
seems to be more natural to start from the usual spin- 
wave picture and expand in fluctuations around an or- 
dered classical gr ound sta te using the Holstein-Primakoff 
transformationPli^IIlEl! Formally, this approach is based 
on the smallness of the inverse spin quantum number 
1/S. However, even for S = 1/2 the result of the 
l/S'-expansion has often been surprisingly accurate, at 
least for quantum antiferromagnets on nonfrustrated lat- 
tices. It is therefore reasonable to expect that the l/S'- 
expansion is a good starting point to describe the mag- 
netically ordered cone state of CS2CUCI4, which is sta- 
ble at sufficiently low temperatures and in a range of 
magnetic fields oriented perpendicular to the layers. It 
turns out, however, that the l/S'-expansion in the cone 
state has been only partially successful because experi- 
mentally observed large quantum renormalizations of the 
exchange energies could not be explained by including 



2 



only the leading \/S corrections to linear spin- wave the- 
ory. In particular, Veillette et a/P calculated the dy- 
namical structure factor in the cone state of CS2CUCI4, 
including the leading l/S corrections to linear spin- wave 
theory. Although they found that the spin-wave interac- 
tions explain, on a qualitative level, many experimentally 
observed features, there was no quantitative agreement 
between theory and experiment. Whether or not infi- 
nite resummations of the 1 / S'-expansion retaining certain 
types of interaction processes would improve the agree- 
ment with experiments remains an open problem. 

The purpose of this paper is to further investigate the 
validity of the spin-wave picture in the cone state of 
CS2CUCI4. Therefore, we present both experimental and 
theoretical results for the magnetic-field dependence of 
the elastic constants and the ultrasonic attenuation rate 
in CS2CUCI4 at low temperatures and in the range of 
magnetic fields where the cone state is stable. Since the 
lattice vibrations are coupled to the spin excitations, the 
magnetic-field dependence of the energy and the damp- 
ing of the phonon excitations provides useful information 
about the magnetic excitations. Indeed, in the vicinity 
of a magnetic phase transition one expects rather strong 
magnetic-field dependence of the phonon degrees of free- 
dom, so that one can map out the phase diagram by 
measuring the elastic constants and the phonon damp- 
ing which determines the ultrasonic attenuation rate.'^S' 

Previous ultrasonic investigations of the longitudinal 
C22-mode on single crystals of CS2CUCI4 by Sytcheva 
et al^ focused on the nature of the different phases 
and on the field- induced quantum critical point. Our 
theoretical approach is complementary to the approach 
adopted by the authors of Ref. [21], who obtained the 
elastic constants from simple phenomenological expres- 
sions for the thermodynamic potentials'^ and relied on 
phenomenological relaxation rates to determine the ul- 
trasonic attenuation rate.'^ In contrast, our derivation 
of the ultrasonic attenuation rate and the elastic con- 
stants is based on a microscopic calculation starting from 
the relevant Heisenberg Hamiltonian. Although the the- 
ory of magnetoelastic interactions is well establishecP^'SlI 
there have been only a few previous investigations of the 
interplay between spins and lattice vibrations in trian- 
gular lattice antiferromagnet^^MUl and other frustrated 
spin systems .fSSl In particular, microscopic investigations 
of the interaction between spin and lattice degrees of free- 
dom in the cone state of CS2CUCI4 have so far not been 
carried out. 

The rest of this paper is organized as follows: Our 
starting point in Sec. |II A| is the assumption that the in- 
terplay between magnetic and lattice degrees of freedom 
in CS2CUCI4 can be described by a spin S = 1/2 quan- 
tum Heisenberg antiferromagnet on a distorted triangu- 
lar lattice with Dzyaloshinskii-Moriya (DM) anisotropy 
in an external magnetic field. The coupling between 
magnons and phonons arises from the fact that the ex- 
change integrals and the DM interactions depend on the 
distances between the actual positions of the ions, which 



in turn depend on the phonon coordinates. Assuming an 
ordere d clas sical ground state (the cone state discussed 
in Sec. II B), we model the low- lying spin excitations in 



terms of suitably defined interacting Holstein-Primakoff 
bosons. By expanding the exchange integrals to second 
order in the phonon coordinates we obtain in Sec. |II C] the 
effective magnon-phonon interaction. In particular, we 
explicitly derive the relevant interaction vertices, which 
exhibit rather complicated momentum dependencies due 
to the spiral spin structure of the cone state. 

In Sec. |III| we then use functional integration tech- 
niques and diagrammatic perturbation theory to calcu- 
late the renormalization of the phonon energies due to 
the coupling to the magnons. Formally, this renormal- 
ization can be described in terms of a momentum- and 
frequency-dependent phonon self-energy, whose real part 
gives the shift in the phonon velocities (which are related 
to the elastic constants), and whose imaginary part gives 
the phonon damping (which is related to the ultrasonic 
attenuation rate). Using the fact that in CS2CUCI4 the 
velocities of acoustic phonons are large compared with 
the magnon velocities, we can derive analytic expressions 
for the magnetic-field-dependent part of the elastic con- 
stants and the ultrasonic attenuation rate. In particular, 
we show that the for small wavevectors k the ultrasonic 
attenuation rate in CS2CUCI4 is proportional to k'^ . 

In Sec. |IV| we present our experimental results for 
the magnetic-field dependence of the elastic constants 
and the ultrasonic attenuation rate in the cone state of 
CS2CUCI4 and compare them with our theoretical pre- 
dictions. In the regime of magnetic fields where our 
perturbative calculation is valid we find good agreement 
between theory and experiment. Finally, in Sec. |V] we 
present our conclusions and discuss some open problems. 
In the appendix we evaluate the contribution of scat- 
tering processes involving intermediate states with one 
phonon and one magnon to the ultrasonic attenuation 
rate; we show that in CS2CUCI4 the smallness of the 
magnon velocities in comparison with the phonon veloc- 
ities guarantees that these processes are small compared 
with the processes considered in Sec. |IIIB[ 



II. MAGNON-PHONON INTERACTIONS IN 
THE CONE STATE OF CsaCuCU 

A. Spin-phonon Hamiltonian 

To model the magnetoelastic properties of the frus- 
trated quantum antiferromagnet CS2CUCI4, we start 
from the following spin-phonon Hamiltonian, 



(2.1) 



where the first term is of the forrrPsEl 

- 2 E [■^'^■^^ • + ■ X ■^^)] - E • 

ij i 

(2.2) 
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Here the Si are spin S = 1/2 operators at positions r^, 
the Jij are the exchange energies connecting spins at po- 
sitions Ti and rj, the antisymmetric vectors Dij model 
the DM interaction, and h is an external magnetic field. 
The second term in Eq. (2.1) describes noninteracting 
acoustic phonons, 



(2.3) 



where ak\ annihilates a phonon with wavevector k and 
polarization A, and the long-wavelength dispersion of the 
phonons is 



WjeA = CA(fc)|A;|. 



(2.4) 



The velocities c\{k) of the phonons depend on their prop- 
agation direction k = The spin-phonon interac- 
tion in Eq. (2.2) arises from the fact that r,j and rj are 



the actual positions of the spins for a given configura- 
tion of lattice distortions {Xi}. Denoting by {Ri] the 
positions of the corresponding Bravais lattice, we have 
Ti = Ri + Xi. The lattice distortions are quantized in 
the usual way, 



X, 
Xk 

Xk\ = 



1 \ " ^tk-Ri y 



kX^kX, 



1 

V^Mujkx 



(akx + a-kx)' 



(2.5a) 
(2.5b) 
(2.5c) 



where the unit vectors Skx define the polarization direc- 
tions of the phonons. Introducing the momentum oper- 
ators conjugate to the Xkx, 



PkX = V 



1 Mcok 



-(flfcA -alkx)^ 



(2.6) 



the pure phonon Hamiltonian (2.3) can be written as 
^P-kxPkx . M 



kX 



2M 



i^kxX-kxX, 



kX-^kX 



(2.7) 



Assuming that positions of the spins deviate only slightly 
from their equilibrium values, we may expand the ex- 
change couplings Jij and the DM vectors Dij in powers 
of the difference vectors Xij = Xi — Xj , 



D, 



J{R^,) + {X,r^r)J{r)l=R 

DiR 



\r=Rij 



(2.8a) 



(2.8b) 



where Rij = Ri — Rj is again a vector of the Bravais 



lattice. Substituting the expansions (2.8a 2.8b) into 
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FIG. 1. (Color online) Triangular lattice with the lattice pa- 
rameters b and c describing the geometry of a single layer 
of CS2CUCI4. The magnetic ions are located at the edges 
of the triangles (black dots) and are mainly coupled via the 
coupling J (solid line) in i-direction and J' along the diag- 
onals (dashed line). The gray dots indicate the position of 
the chlorine atoms mediating the exchange interaction, which 
depends in a simple model on the three distances h, I2, h and 
on the two bonding angles ai and 02 • For simplicity, we lump 
these dependencies into the variation of the bond length. In 
the same manner, also the dependence of the coupling J' on 
the angle (p is neglected. 



Eq. (2.2) we obtain an expansion of our spin-phonon 



Hamiltonian in powers of the phonon operators. 



H- 



pho 



H, 



Ipho 



H. 



2pho 



(2.9) 



where the pure spin part i/spin is obtained by assuming 
that the spins are located at the sites Ri of the Bra- 
vais lattice, while the spin-phonon interactions H^p^^ and 

^s^pin° describe the coupling of one and two powers of the 
phonon operators Xij to the spin system. 

The pure spin part -ffgpin of our spin-phonon Hamil- 
tonian (2.9) which is formally identical with Eq. (12. 21) 



except that now the spins are located at the sites Ri of 
the Bravais lattice. From the values compiled in Tab.|l]it 
is clear that for the layered material CS2CU CI4 th e cou- 
plings between neighboring layers is very smalP^ZI go that 
we may ignore the interlayer coupling J" and focus on 
a single layer. The spins are then located at the sites of 
a distorted triangular lattice characterized by the lattice 
parameters b and c, as shown in Fig. [T] Within the planes 
of the lat tice, we assume an anisotropic nearest-neighbor 
exchange ,1^^^ 



Jij — J i^Rji R'j^ 



where 



J if R, - R^ = ±{Si + 62) 
J' if R, - Rj = ±di or ±^2 ' 

(2.10) 

(2.11a) 
(2.11b) 



b c 

01 = -e^ + -By, 

02 = -j^x - -je.y, 



so that 82= be,j. (see Fig. [ij). The DM vectors Dij 
are assumed to point in the z-direction and connect sites 



Parameter 


Value (meV) 


J 


0.374(5) 


J' 


0.128(5) 


J" 


0.017(2) 


D 


0.020(2) 



27r/g, 



TABLE I. Accepted valuefP^l 

of the in-plane exchange in- 
teractions J and J', the inter-plane interaction J", and the 
Dzyaloshinskii-Moriya interaction D in CS2CUCI4. All com- 
parisons between theory and experiment presented in this pa- 
per are based on these values. Note that the largest exchange 
coupling J = 0.374 meV corresponds to a temperature of 4.34 
Kelvin. 



in the directions ±<5i and ±^2, 

= D{R,-Rj)e, 



(2.12) 



with Dij = —Dji given by 



D, 



D{R, -Rj)^ ±D if R, - Rj ^ ±5i or ±82. 

(2.13) 

For CS2CUCI4, experimental estimates for J, J', D, and 
the interlayer exchange coupling J" are summarized in 
Tab. |l] Throughout this paper, we assume that the ex- 
ternal magnetic field is perpendicular to the plane of the 
lattice, h — he^. In the experimentally relevant regime 
J' < 2J the spin system has then a unique classical 
ground state, as discussed in the following section. If the 
magnetic field lies in the plane of the lattice the ground 
state is more complicated and cannot be described within 
the framework of the spin-wave expansion used in this 
paper.E^ 

B. Spin- wave expansion 

1. General procedure 

To set up the spin- wave expansion, we should first iden- 
tify the ground-state spin configuration in the classical 
limit, where the spin operators are replaced by three- 
component classical vectors of length S, i.e.. Si — > Srhi, 
where rhi are unit vectors which point in the direction of 
the local magnetization.!^ In the experimentally relevant 
regime J' < 2 J the classical ground state of our model is 
a spiral in the plane of the lattice, which is tilted towards 
the direction of the magnetic field. In this so-called "cone 
state" the magnetization points locally in the directiorP 

rhi = ce[cos(Q • Ri)e,^ + sm{Q ■ Ri)ey] + se e^, (2.14) 

where we have introduced the abbreviation 

sin6' , ce = cosd. (2.15) 



This state is characterized by the opening angle 9 of 
the cone and the wavevector Q of the spiral, as shown 
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FIG. 2. (Color online) Projection of the graph of the "cone 
state", which is the classical ground state of the anisotropic 
triangular lattice antiferromagnet for J' < 2J and magnetic 
field perpendicular to the lattice plane. Note that the wave- 
length 2n/Qx of the spiral has the indicated length and points 
along the i-axis in the direction of the arrow. 



in Fig. [2] To bosonize our spin Hamiltonian using the 
Holstein-Primakoff transformation, we follow the general 
procedure outlined in Refs. 28 and [29] and complement 
the unit vectors rhi by two orthogonal unit vector e|^-* 

and ep'' such that {ef^\ e'f \ rhi} form a right-handed 
triad of orthogonal unit vectors for each site Ri. The 
transverse basis vectors el^-* and e^^'' are defined only 
up to a local U{1) gauge transformation.f^Sl A convenient 
choice is 

e-^-* — sm{Q ■ Ri)ex — cos(Q • Ri)ey, (2-16) 

(2) 

el = se[cos(Q • Ri)e^ + sm{Q ■ Rt)ey] - cg e^. 

(2.17) 

Next we introduce the corresponding spherical basis vec- 
tors 



,(2) 



^i = e] ' +ipel ' , p = ±, 
and expand the spin operators in this basis, 

s, = s\7h, + sl, 

with the transverse part given by 



(2.18) 
(2.19) 

(2.20) 



The spin components are then bosonized using the 
Holstein-Primakoff transformation. 



S - Hi 



S+ ^V2Sjl-^h, 



sr ^V2SbU i- 



rii 
2S' 



(2.21a) 
(2.21b) 

(2.21c) 
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where rij = blbi and bi are canonical boson operators. 
The bosonized spin part of our Hamiltonian can then be 
written as 



where 
Ho 



spin =Ho+ i?2|| + +H^+ i/ll_L, 



(2.22) 



(2.23) 



S 



H2\\ =-^J2 4("* + '^j)+^f^- (2-24) 



Hm = 



9 



(2.25) 



2 z_^^u'^i'^i' 
-gEE^''^^'^/''' (2-26) 

ij pp' 

= E ^'^ • [(-^'J - AjX)m,(5 - n,) - S,jh] 

ij 

= -Y, ■ - D., X )rh,n,] , (2.27) 



where in the last line we have used the fact that the spin 
configuration in the classical ground state satisfie^^U 



m," X 



= 0. 



(2.28) 



For simplicity we have introduced the following effective 
couplings between longitudinal and transverse spin fluc- 
tuations, 

jfj = Jijifii ■ rhj + ■ [rhi x rhj), (2.29a) 

J^'=J,,(er.ef) + A,-(ef xef). (2.29b) 

Defining (pi = Q ■ Ri and cpij = tpi — tpj = Q ■ {Ri ~ Rj), 
we obtain with our choice of basis 

"^'i = "^ij [^e + cos ipij] - D^jcj sin ip^j , (2.30a) 

= Jij[c^g + (1 + sl) COS(^ij - 2isg siliLpij] 

— Dij [(1 + si) sin ipij + 2isg cos ipij] , (2.30b) 

= -cg [ (1 - cos ipij) + Dij sin ipij] . (2.30c) 

It is convenient to introduce the Fourier transforms of 
the exchange couplings and the DM interaction: 



R 



= 2Jcos(fc^fo) + 4J'cos(fc^6/2)cos(fcj^c/2), (2.31) 
Dk = J2e-'''''DiR) 



R 



-UD sin(fcj;6/2) cos(fcyc/2). 



(2.32) 



We also define the combination 

Jk=Jk-iDk, (2.33) 

which is real because the Fourier transform of the DM 
coupling Dk is purely imaginary. The Fourier transforms 
of the effective couplings (2.30a -2.30c I can then be writ- 
ten as 



Jh. — ^eJk 



jD , jD 
2 "^Q+k + "JQ-k 



cpk + {l + s 



2 

•J r 



2-, "Q+k -JQ-k 



+ [Jq+k - Jq-k] 



■^k 



— ^k — 



Jk — 



jD I jD 
•jQ+k ^ '^Q-k 



(2.34a) 

(2.34b) 
(2.34c) 



2. Classical ground state: Cone state 

To fix the parameters 6 and Q which characterize the 

into our 



classical ground state, we substitute Eq. (2.30a 



expression (2.231 for the classical ground-state energy 



and minimize with respect to 9 and Q. The classical 
ground-state energy can then be written as 

Ho = N^jl^^-NShsg 

Minimizing this with respect to 9 we obtain 

se=sii^9 = h/hc, (2.36) 

where the critical magnetic field is given by 

K = S{J^ - J^) = SiJo -Jq + iDq). (2.37) 

The wavevector of the spiral is obtained by minimiz- 
ing Hq in Eq. (2.35) with respect to Q, which amounts 
to finding the minimum of the real function Jq. The 
wavevector of the spiral is thus determined: 



Jq = Vfc {Jk - iDk)k=Q = 0- 



(2.38) 



Anticipating that this condition leads to a spiral along 
the a;-axis, Q = Qx^x, and using the above expressions 



for Jfc and Dk, it is easy to show that Eq. (2.38) reduces 



to the following transcendental equation for x = Qxb: 



cos 



(I) 



J' D (X 

cot — 

2J 2J \2 



(2.39) 



Note that for Z) = and isotropic exchange couplings 
(J' = J) this condition reduces to cos(a:/2) = —1/2, im- 
plying X — Qxb — 47r/3, which describes the usual 120° 
ground state of an isotropic triangular lattice antiferro- 
magnet. 
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3. Magnon dispersion 

To calculate the magnon spectrum to leading order in 
the 1/5 expansion, we may approximate ~ V2Sbi 
and 5" w \/2Sb\, so that the transverse part of our 
spin-wave Hamiltonian is approximated by H±^ w -ff2_L 
with 



j 



(2.40) 



The magnon spectrum can then be obtained by diagonal- 
izing the quadratic boson Hamiltonian H2 — i?2|| +^^2±- 
Combining Eqs. ( |2.24[ ) and ( |2.40[ ) we obtain 



+ J^.^bjb] + Jr^hb,] ]+hseY, (2-41) 

i 

Introducing the Fourier transform of the boson operators 
via 



we obtain in momentum space, 
^[Akb'^bk+^ 

k 

where Ak = + , with 



H2 = + ^ W-k + b-kbk]}. (2.43) 



A1. 



-Bk — S 



J, 



Q 



Bk = — — Cfl 



jD _ jD 



Jk — 



jD , jD 
•^Q+k ^ -^Q-k 



(2.44a) 

(2.44b) 
(2.44c) 



These coefficients are real and have the symmetries 
A'^f, = ±A^ and Bk = B_k- Note that 



S 

Bo = --^cg[Jo - Jq] = --^hc, 
Ao = -^K, 



(2.45a) 
(2.45b) 



which implies that the magnon dispersion is gapless at 
fc = 0, as required by the C/(l)-symmetry of the Hamil- 
tonian. Note also that for small fc the linear coefficient 
in the Taylor expansion o f the antis ymme tric coefficient 
A^ vanishes due to Eqs. (2.33) and (2.38), so that 



Al 



h = O.l/i. 
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FIG. 3. (Color online) Graph of the magnon dispersion Ek of 
the anisotropic triangular lattice antiferromagnet CS2CUCI4 
with J'/ J = 0.34, Dzyaloshinskii-Moriya anisotropy D/J = 
0.054, and for a magnetic field h = O.lftc; see Eq. (2.50 1. Note 



that the f/(l)-symmetry of the Hamiltonian guarantees that 
for fe = the magnon dispersion is gapless.'^ 



To obtain the magnon dispersion, we diagonalize the 
Hamiltonian (2.43) using the following canonical (Bogoli- 



(2.42) ubov) transformation. 



bk 

6L 



Uk -Wfe 
-"^fe Uk 



Pk 
Plk 



(2.47) 



where 



Uk 



2€k 



Vk 



Bk A+- Ek 



and 



\Bk\ 



^k = ^{AlY-Bl. 



2efc 



(2.48) 



(2.49) 



The energy = e„fc is even under fe — > — fc; the full 
magnon dispersion is 



Ek — ^k + A. 



{At? 



Bl + A,, 



(2.50) 



which does not have any definite symmetry with respect 
to fc — > — fe. The diagonalized form of the magnon Hamil- 
tonian reads 



E 



Ek/3i/3k + — 



At 



(2.51) 



A graph of the magnon dispersion _Efc is sh own in Fig. [3] 
Using the fact that according to Eq. (2.46) the term A^ 
can be neglected for small fe, the leading long- wavelength 
limit of the magnon dispersion is 



(2.46) Ek = Jv^,kl + vlkl + Oik')^vik)\k\+Oik'), (2.52) 
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FIG. 4. (Color online) Plot of the linear spin-wave result 
for the magnon velocities of CS2CUCI4 in the two principal 
directions of the lattice as a function of the magnetic field. 
The relevant values of J, J', and D are given in Tab.jl] It turns 
out that for an arbitrary magnetic field the magnon velocities 
are small compared with the velocities cx of acoustic phonons 
in CsaCuCU. 



with direction-dependent magnon velocity 
v{k) = 



vlkl, 



(2.53) 



where fc = fc/|fc|. In Fig. [i] we show a graph of the linear 
spin-wave result of the two principal magnon velocities 
Vx and function of the external magnetic field for 

CS2CUCI4, where J, J', and D have the values given in 
Tab.m 



4- Interactions between magnons 

The terms in the transverse Hamiltonian H±^ involving 
more than two boson operators and the contribution s H^\ \ 
and i?||_L to the magnon Hamiltonian given in Eqs. ( 2.25 ) 



and (2.27) describe various interaction processes between 
magnons. As compared with the S'-dependence of the 
free magnon dispersion on S, the interaction contains 
higher orders of the small parameter 1/S, so that at first 
sight it seems that for large S we can simply neglect inter- 
action effects. It turns out, however, that for the proper 
calculation of the long- wavelength limit of the ultrasonic 
attenuation rate the leading contribution H3 to the mix- 
ing term ffyj. involving three magnon operators cannot 
be neglected. This part of the magnon Hamiltonian can 
be written as 



-Cg- 



2i 



E 



,b]b,h 



(2.54) 



where we have defined 

Kij = S0 [Jijil - cos (fiij) + Dij sin ipij] 
+ i [Jij sin (fiij + Dij cos (pij] . 

In momentum space Ht, can be written as 
1 

~ fcifc2fc3 



(2.55) 



- E ^fcl+fc2+fc3:0 



^Tt^'\k,M-M)hUb\b^, 



l^Tt''\k^-M,k^)bl^bk,bk, 



fci 



(2.56) 



with the vertices given by 



rf^'^(fei,fe2;fc3) 



-C0- 



2i 
2S 



[K.k,+K_k,], (2.57a) 



(2.57b) 



.3 ^. 
where is the Fourier transform of the function Kij 



defined in Eq. (jS^S]) 

Kk = Sg Jk - 



jD , jD 



jD „ jD 



For later reference we note that 



Ka = Sg{jQ - Jq) = Sg-^. 



(2.58) 



(2.59) 



C. Magnon-phonon interactions 

1. Modeling the derivatives of the exchange couplings in 
CS2 CuCk 



Substituting the gradient expansions (2.8a I and (2.8b) 
for the exchange and DM coupling s int o the spin-phonon 
Hamiltonian defined in Eq. (2.2 1, we obtain an ex- 



spm 



pansion of H^^-^^ in powers of the phonon operators of 



the form (2.9). Since the exchange interactions are large 
compared with the DM interactions, we shall neglect the 
dependence of the DM couplings Dij in the phonon co- 
ordinates. The n-phonon part of our spin-phonon Hamil- 
tonian can be written as 



spin 



^^U^;>S,-S,, (2.60) 



where the coupling functions involving n = 1 and n = 2 
phonons are 



(1) 



(1) 



(2) 



(2.61) 



(2.62) 
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with the vector jj^p and the tensor 3^^^ defined by 



J(r) 



(2.63) 
(2.64) 



r=Rij 



The Fourier transforms of functions (2.61) and(2.62) are 



1 



r(i) 



J 



fc' ^ 



(2.65) 



2iV 



where 



-(2) _ }_\^ -ik-Ri-ik'Rjjj{2) 



jW^^e-^'^-«j(i)(i?), 



(2.66) 
(2.67) 



(2) 

and similarly for . Note that the Fourier transform 
j'jp of the gradient has the symmetries 

which implies that j'^^ is purely imaginary and vanishes 
for fc = 0. In contrast, the second derivative tensor is an 
even function of k. 



(2.68) 



r(2) 



r(2) 



(2.69) 



and thus has a finite limit at zero momentum. In general, 
in-plane lattice deformations in CS2CUCI4 induce changes 
in both couplings J and J' as indicated in Fig. [l] In 
principle one should use four independent parameters to 
describe the changes of J and J' in response to the dis- 
placement of the lattice sites in two independent in-plane 
directions. To reduce the parameter space, we note that 
usually the exchange couplings originate from an elec- 
tronic exchange path where orbitals of neighboring atoms 
overlap. For CS2CUCI4 the exchange is mediated by the 
chlorine atoms located in between the copper atoms, as 
shown in Fig.[T] In a rough approximation where only the 
coupling of a few atoms is considered, the exchange cou- 
plings do not change to first order in the displacements 
when one of the copper atoms moves up or down. There- 
fore one can neglect the variation of J with respect to 
displacements in the y-direction and assume that J only 
changes appreciably if the atoms move in the x-direction. 
By a similar argument one can also neglect the depen- 
dence of the couphng J'(r) = J'(r, ip) on the polar angle 
if, setting 



dJ' 
dx 



cos (po 



djr_ 

dr 



dJ' 
dy 



sm ipo - 



dJ' 



(2.70) 



where sintpo = c/{2d). Here d = vP + c?/2 is the rele- 
vant bond length, and dJ' /dr is the derivative of the cou- 
pling J' with respect to the distance between the spins. 
With these approximations we obtain the following ex- 



pression for the gradient 
in CS2CUCI4: 



(1) 



of the exchange couplings 



r(i) 



4i 



^ dJ . /, ,N 
.dJ' 



dr 



Ex cos ifo sin {kxh/2) cos {kyc/2) 



+ eySunp()COs{kj.h/2)sv[\{kyc/2) . (2.71) 



Similarly, we obtain for the matrix elements of the second 



(2) (2) 
derivative tensor ^ = ^ , 



[jif^]^^=2cos(A:,6) 



9V 



9V 



-I- 4 cos^ (yso cos{kxb/2) cos{kyc/2) -^-^ , (2.72a) 



"^k^] 12 ~ ''^^ 'f'^ ^^"^ sin{kxh/2) sin{kyc/2) 



Qj.2 ' 

(2.72b) 



^2 J/ 

[Jfc ^] 22 = 4 sin^ ipo cos{kxb/2) cos{kyc/2) . (2.72c) 



Due to the [/(l)-symmetry of the Hamiltonian, the vec- 
tor J^^' and the tensor jj,^' are only fixed up to a phase 

that depends on the choice of triad {e'"^\ ep\ mi}; see 
Ref. [28] for a discussion of this point. Note that the ge- 
ometry shown in Fig. [l] suggests that the couplings will 
become stronger if the bond length is enlarged; at first 
sight this counterintuitive behavior originates from the 
fact that the bond angles between the chlorine atoms be- 
come larger with increasing bond length, which usually 
results in a stronger coupling.^Sl Later we will see that 
only the squares of the derivatives enter in the final re- 
sult for the velocity shifts and the ultrasonic attenuation 
rates, so that one cannot directly determine the sign of 
the change of the coupling strengths from the ultrasound 
measurements . 



2. Phonon renormalization in the classical spin background 



Expanding the spin operators in Eq. (2.60) in powers 



of Holstein-Primakoff bosons, we obtain a double expan- 
sion of the spin-phonon Hamiltonian in powers of phonon 
operators and magnon operators. For our purpose, it is 
sufficient to retain only terms up to two phonon opera- 



tors. After expressing the spin operators in Eq. (2.601 in 



terms of Holstein-Primakoff bosons, to leading order for 
large S', we generate also pure phonon contributions of 
the form 



H, 



npho 



(2.73) 
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which describe motion of n phonons in the classical spin 
background. In particular, the one-phonon term Hq^^° 
can be written as 



Ipho 

~ 



2 



U' 



(1) 

Q.Q 



(1) 

Q. Q 



(2.74) 



(1) 
Q.Q 



But from Eqs. (|2.65|) and (|2.68|) we see that C/^ J 



0, implying Hq^^^° = 0. On the other hand, 



the corresponding two-phonon part 

r.2 



H, 



2pho 




2 



U 



(2) 

Q,Q 



U. 



(2) 

Q-Q 



(2.75) 



is finite and yields an important contribution to the mag- 
netic field dependence of the elastic constants. Inserting 
the explicit expression for the Fourier transform U^, j^, of 



(2) i I, 

the second derivative couplings U-,' [see Eqs. (2.62) and 



( 2.66[)] an d using the symmetry Eq. (2.691 we see that 
Eq. ( |2.75 l can be written as 



„2pho ^ 



kX 



E[j''°(fc, A)X_feA-'^feA 



(2.76) 



where the zeroth order contribution to the phonon self- 
energy (in an expansion in powers of the magnon opera- 
tors) is given by 



^pho 



(fc,A) 



^ t 



Here we have defined 



t(2+) _ t(2) 



r(2) 



,(2) 
'k+Q 



r(2) 



r(2) 
'k Q 



^e'^Q,k 



e-kx ■ 
(2.77) 



k,Q 2 

For later convenience we also introduce the notation 



(2.78a) 



,(2-) 
k,Q 



r(2) 
'k+Q 



,(2) 
'fc-Q 



(2.78b) 



Using 

ef ■ rhj = C0 sin ip^j =F icese [1 - cos ipij] , (2.81) 
and Fourier transforming to momentum space, we obtain 



H 



Ipho _ ^ pX6 
k 



with the hybridization vertex given by 



Here we have defined 



- J 



'k.Q 



r(l) 


•^k+Q 




'^k-Q 


Q 




2 




r(l) 


'^k+Q 


+ 


^k Q 


fe 




2 





(2.82) 

(2.83) 

(2.84a) 
(2.84b) 



For fixed Q the functions J^q' have the symmetries 



r(l±) 
-k,Q 



(1±) 
k,Q ■ 



(2.85) 



The contribution (2.82) to the Hamiltonian describes the 
hybridization between magnon and phonon modes. Note 
that for finite Q the relevant coupling vanishes lin- 
early in fc for fc — > 0. As it will be discussed in Sec. |III A{ 
the hybridization term -ffi^ag gives rise to an additional 
renormalization of the phonon velocities which has the 
same order of magnitude as the renormalization due to 
the contribution SQ^°(fc, A) aris ing from the classical spin 
background given in Eq. (2.77). 

Note that the vertex F^'' does not have a definite sym- 
metry under A; — > — fe. To obtain a more transparent 
classification of the vertices, it is useful to express the 
Holstein-Primakoff bosons bk in terms of Hermitian op- 
erators ^k = and Hfc = h1^ as foUowsPSSIl 



(2) 

Keeping in mind that is an even function of fc, we 

(24-) (2 — ) 

see that q is even while q is odd under fc — > — fc, 



t(2±) _ I t(2±) 



(2.79) 



3. One-phonon one-magnon hybridization 

The part hI^^"^ of our spin-phonon Hamiltonian in- 
volving one phonon and one magnon operator is given 

by 



rrlpho _ 
linag 



{2S) 



3/2 



b,{e^ -rhj) . 

(2.80) 



2 ,/2Afc 



n-fc. 



where the energy is given by 



Ak = A+-Bk = A+ 



\B, 



(2.86) 
(2.87) 

(2.88) 



Then our quadratic magnon Hamiltonian (2.43) can be 
written as 

H2 = \ ^{n^fcHfe + e^$_fc$fe 
fe 



tA^ ($_fcHfc-$fcH_fc)-A+|, 



(2.89) 
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while the magnon-phonon hybridization (2.82) reads 



where the couphng functions are 



k 



i k 



(X^fcHfe-XfeH^fc)}, (2.90) 



where the vertex F^* describing the couphng of the 
phonon coordinates to the generahzed magnon coordi- 
nates $fc is an even function of fc, while the correspond- 
ing vertex T^'^ which couples the phonon coordinates to 
the magnon momenta is odd in k: 



^ k 



^ k 



^k r-r^Xh 



~k ■ 



(2.91a) 
(2.91b) 



In principle, the complete quadratic part H2 



HP 



ho 



H, 



2pho 



^imag of the magnon-phonon Hamil- 
tonian can be diagonalized by means of a canonical 
transformation.'^'^ But due to the absence of inversion 
symmetry of the magnon dispersion the explicit con- 
struction of this transformation seems to be a rather com- 
plicated mathematical problem which we have not at- 
tempted to solve. Fortunately, the renormalized phonon 
spectrum can be obtained in a straightforward manner 
at this level of approximation by integrating in the path 
integral formulation of the theory over the magnon field. 
As will be shown in Sec. Ill A[ the phonon spectrum can 
then be directly obtained from the inverse phonon propa- 
gator appearing in the resulting Gaussian effective action 
for the phonon field. 



4. One-phonon two-magnon interaction 

The part of the spin-phonon Hamiltonian involving one 
phonon and two magnon operators is 



+ l[{et-ej)blb, + {e--el)b,b] 
+ (e+ . e+)blb] + (e- ■ ej)b,b,]}. (2.92) 



Using Eqs. (2.16 2.18) to calculate the relevant scalar 
products of the basis vectors, we obtain after Fourier 
transformation to momentum space. 



H. 



Ipho 
2mae 



^k.k' ■ Xk-k'blbk' 



-y 

^ (rfc,fc' • Xk+k'blbl, + Tfk' ■ X-k-k'bkbk' 



(2.93) 



■ k,k' 



2c; 



J, 



- ( r ■ 

k-k',Q y'^k,Q 



_ t(1-)> 



Tb*b^ _ rbb _ 
^ k.k' — — A t 1,' — 



s 



k,k' 



t(I-) 



Since the Fourier transform J^"^ of the derivative of the 
exchange coupling is purely imaginary, this is also true for 



, j(i-) 

k'Q 



(2.94) 
(2.95) 



Note 



the above vertex functions T^'^j.,, rf].! and T^'j^, 

also that the normal vertex ^, is antisymmetric while 

the anomalous vertices F^*^^, and F^''^., are symmetric 
under k k' . Keeping in mind that these vertices are 
purely imaginary, one easily verifies the following rela- 
tions, 



(2.96a) 
(2.96b) 
(2.96c) 



Of particular interest is the leading behavior of the above 
magnon-phonon vertices for small k and fc'. 



j^b^b 
^ k,k' 


Y^bH 
— k',k — 


i^ll'%r 




^ k,k' 


— ^ k'.k — 


j^bb 

^ k',k — 


iri\kr 


j^bb 
^ k.k' 


j^bb 

— ^ k',k — 


^ k',k — 





k.k' 



s 



2 ^i2sl-cl)[{k-k')-Vci][j^^^ 
se[{k-VQr^ik'-VQf]j^^'^^ 



-J. 



1)1 



S 



k.k' 



cl[{k + k')-VQ][j'^ 



(1) 



Q=0 



Q J 
0(fc3), 

(2.97) 
+ 0(fe3). 
(2.98) 



It turns out that -ff2mag El'l- (2.93) is not the only 
part of the Hamiltonian describing tlie coupling of a sin- 
gle phonon to two magnon operators. The reason is 
that due to the magnon-phonon hybridization discussed 
in Sec. Ill C 31 the Holstein-Primakoff bosons are a linear 
combination of true magnon operators and phonon oper- 
ators. In principle, the proper linear combination can be 
obtained by diagonalizing the quadratic magnon-phonon 
Hamiltonian H2 + HP^'° + H^'^'^° + hI^°^ by means of a 
canonical transformation. Although the construction of 
such a transformation is in principle possible,'^ for our 
purposes it is fortunately not necessary to explicitly solve 
this complicated algebraic problem. The reason is that at 
long wavelengths and to the order in the small parameter 
1/S consistent with Eq. (2.92 1, the proper magnon oper- 
ators &fc can be determined from the requirement that in 
the quadratic part of the magnon-phonon Hamiltonian 
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J?2 + HP^° + H^^^° + hI^°^ there should be no terms 
describing the couphng of the phonon coordinates Xk to 
the generahzed magnon momenta 11^, as in the second 
hne of Eq. (2.90). It is then easy to show that the proper 



magnon operators are related to the original Holstein- 
Primakoff magnons bk via the phonon-dependent shift 
transformation 



where 



Afc — 



Xk 



-.xn 



(2.99) 



(2.100) 



wavelengths, 

STii,=-2Ssl[{k-k')-VQ][j^^^ 

<t'=^.s^[(fc + fc')-VQ][4^' 
+ 0{k^). 



Q=0 



J. 



(1) 



Q=0 



J. 



(1) 



(2.104) 



(2.105) 



Combining these expressions with the long-wavelength 
limits of the corresponding bare vertices given in 
Eqs. (2.97 2.98), we find that the proper one-phonon 



two-magnon vertices are in the long-wavelength limit 
given by 



depends on the antisymmetric vertex defined in 

Eq. (2.91b I. If we now substitute the transformation 



(2.99) into the part Hy, of the pure magnon Hamiltonian 



involving three powers of the Holstein-Primakoff bosons 



given in Eq. (2.56), we generate (among other terms) an 
additional contribution to the one-phonon two-magnon 
interaction, which involves the same power of S as i?2mag 



^ k,k' 

s 



[2sl + cl)[[k-k')-yQ\[j, 



'q 



se [{k ■ Vq)2 - (fe' • Vq)2] J« l + 0(fc3), (2.106) 



S , 



in Eq. (2.92| and therefore should be taken into account Tfc,!' — "^(^Sg -I- Cg) [(fe + fc') • Vq] [ Jq'' 



on equal footing with -ff2mag- Fortunately, this contri- 
bution can be absorbed via a simple redefinition of the 



vertices in Eq. (2.93) 



^ k,k' — fc,fc' 

^ k,k' — ^ k,k' 

^-^bb T^bb 

^ k,k' — fc,fc' 



' ^^k.k'^ 
' ^^k',k' J 



(2.101a) 
(2.101b) 
(2.101c) 



where the correction terms due to magnon-phonon hy- 
bridization are 



2 

0(fc3). 



Q=0 



(2.107) 



Note that the coefficients of the terms linear in the mo- 
menta differ only by a minus sign, which will turn out to 
be essential to obtain the correct long-wavelength limit 
of the ultrasonic attenuation rate. 

At this point it is convenient to express the effective 
one-phonon two-magnon Hamiltonian -ff2mag terms of 
the Bogoliubov quasiparticle operators (3^ and defined 
in Eq. (|2.47|). We obtain 



Srtl ^Vi'^'ik - k' ,-k;k')Xk--k 



Eg ''''(— fc; fc', k — k')\k-k' 



6r 



k.k' 
-.66 



-^6+6+6 



(— fc, ~k'; k + k')\k+k' 



6T''lk'=n'\-k-k';k,k')\k+k' 



(2.102a) 
(2.102b) 
(2.102c) 



The effective one-phonon two-magnon Hamilt onian 



^2mag can then be obtained from i?2mag '^^ Eq. (2.93) 



Ipho 



by substituting the bare vertices by the shifted ones. 



7-16*6 , -r6t6 T-v6t6+ 
^ k,k' ~^ ^ k,k'' ^ k,k' 



rfj Und nV^ fg^,. Using 



the fact that according to Eqs. ( 2.57a 2.59 1 at long wave- 



lengths the three-magnon vertices can be approximated 
by 



rf '^(fci, fe2; fea) = -rf'^i-k,, -kr, -fci) 



-CeSg- 



25 he 



0{k^), (2.103) 



we obtain for the shifts of the one-phonon two-magnon 
vertices due to magnon-phonon hybridization at long 



H. 



Ipho 



E 

kk' 



^k.k' ■ Xk-k'Pil3k' 



21 {^kl' ■ ^k+k'l^lPl, + ■ X^k-k'(3kPk' 



where the vertices are given by 



k,k' 



k,k' 



k,k' 



f,6+6 ■i^6t6 

= UkUk'Tk^k' - VkVk'i--k-k' 

UkVk'tt'lk'-VkUk'f\k', 

-ukUk'ti'k' + vkVk'r''X-k' 

rb^b I fi6t6 

UkVk'Tk^^k' + VkUk'T ^k.k' ^ 



^ ^'^kUk'T'lk' + VkVk'T^X-k' 
UkVk.ttlk' - VkUk'f'-Xk, 



(2.108) 



(2.109a) 
(2.109b) 
(2.109c) 



As will be shown in Sec. |IIIB[ at zero temperature the 
ultrasonic attenuation rate is determined by the anoma- 
lous vertex f ^ ^, , whose long-wavelength limit is explic- 
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itly given by 



^ k.k' 



24 

2 



1 



[(fc + fe')-VQ][j, 



2 sj^h^k' 



[(fe-VQ)2~(fc'.VQ)2]j^^ 



Q=0 

(1) L I /o/i-s 



(2.110) 



III. PHONON SELF-ENERGY DUE TO 
MAGNON-PHONON INTERACTIONS 

A. Elastic constants 

The elastic constants are directly related to the veloc- 
ities of the acoustic phonons which can be determined 
experimentally with high accuracy!^ Although the spin- 
phonon coupling is expected to be small, its influence 
on the phonon properties is visible in the magnetic- 
field dependence of the elastic constants. The leading 
contributions to the shift in the elastic constants is al- 
ready contained in the quadratic magnon-phonon Hamil- 



2pho 



tonian E-^ + -f H'^ 

energy Yj^°{k, A) given in Eq. ( |2.77 



Ipho 
I mag • 



The phonon self- 
) which is due to 

the coupling of the phonons to the classical spin back- 
ground renormalizes the phonon frequencies according to 
-Ti^°{k, A). The self-energy correction leads 



■"kx 



to the following relative shift of the phonon velocities, 



lini 
\k\^o 



1 



(3.1) 



An additional renormalization of the phonon veloci- 
ties arises from the magnon-phonon hybridization in 
Eq. (2.82). To calculate this contribution, wc note that 



the phonons couple to the magnons only via the com- 
bination Xk = X^A -^kx^kx SO that it is convenient to 
describe the phonon dynamics within an effective La- 
grangian which is obtained by integrating over the canon- 
ical phonon momenta Pkx in the phase space functional 
integral representation of the theory. ■^^ To leading order 
for large S we may truncate the effective Euclidean action 
at the quadratic order in the fluctuations, 

S[X, p, 13] « 52pho[^] ^ 52mag[/3, /?] + ^l^a^^- ^> 1^1 

(3.2) 

where the Gaussian actions S'^p'^°[X] and S'2mag[/3, /?] de- 
scribe noninteracting phonons and magnons, 

52pho[x] ^^Y. ^^[^' + + (fc, A)] 



KX 



X X_ 



kxXkx , 



2inag 



K 



K 



(3.3) 
(3.4) 



Here T is the temperature and K = (iw, k) is a collective 
label containing bosonic Matsubara frequencies iuj and 
wavevectors k. The real field Xkx represents the Fourier 
components of the phonon operator Xk\, while the com- 
plex boson field 13k represents the Fourier components of 
the magnon operator (3^. From Eqs. (2.47 2.82) we see 



that the magnon-phonon hybridization is represented by 
the Euclidean action 



/3] - ^ E " • {X-kPk + Xk 

K 



(3.5) 



where the hybridization vertex is the following linear 
combination of the corresponding hybridization vertex 
in the Holstcin-Primakoff basis given in Eq. (2.831, 



^ k 



^k^k 



(3.6) 



The Gaussian integral over the magnon field is now easily 
carried out, and we obtain for the effective phonon action 
in Gaussian approximation. 



KXX' 



<5A,A'M[c.2+c.2^ + EP''°(fc,A)] 

X-kxXkx' ■ (3.7) 



iix) — Ek 

For simplicity, we neglect in the sum over the phonon 
modes the optical phonons and off-diagonal terms A 7^ A'. 
In this approximation the magnon-phonon hybridization 
generates the following phonon self-energy: 



^f°{K, A) 



fifcA 



Miiuj - Ek) 



(3.8) 



Combining this with the classical self-energy Yi'^°{k,\) 
given in Eq. (2.771, we conclude that the renormalized 



phonon dispersion uj^x is determined by the real positive 
root of the cubic equation 



wIa 



^feA 



^feA 



^plio 



(fe,A) 



T.f\k,X) 



^pho / ■ 

[lU! 



Skxl 



u}kx,k,X) 

|2 



M{{bkx-Ek)' 



(3.9) 



Recall that according to Eq. (2.77) the classical self 



energy I]Q^°(fe, A) is proportional to . Keeping in mind 



that according to Eq. (|2.83|) the hybridization vertex 

is of order S*^/^ and the denominator in Eq. (3.8 ) contains 



one power of E^ <x S, we see that both self-energy con- 
tributions on the right-hand side of Eq. (3.9 1 have the 



same order of magnitude in a formal 1/5-expansion. 

Because in the experimentally relevant regime the 
magnon velocity is small compared with the phonon ve- 
locity, the shift in the phonon velocities due to magnon- 
phonon hybridization can be approximated by 



(Aca)i 
Ca 



lim - — „ , ^ 
|fc|-i.o 2Muj 



GfcA 



(3.10) 



fcA 
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Note that this contribution is always positive, i.e. the 
couphng to the magnons enhances the phonon veloci- 
ties, which is a simple consequence of the fact that for 
CS2CUCI4 the magnon velocities are small compared with 
the phonon velocities. To explicitly take the limit in 
Eq. (3.10) we need th e sma ll-mo men tum limit of the ver- 
tejiTf, . Using Eqs. (2.83) and (3.6) we obtain to leading 
order 



-(25)3/2|fe|3/2. 



(3.11) 



where the dimensionless vector F-^^{k) is given by 

i^^^(fc) = sef^^{k) - clf^^ik k) , (3.12) 
and we introduced the auxiliary vector functions 



/f^(fc) = -(fc.VQ) 



1 



Q=0 



J, 



(1) 

Q 



f^'\k,k') = — ^(fe . VQ)(fc' • Vq)J^^ 

2v[kj 



(1) 



(3.13a) 
(3.13b) 



Here we have approximated the magnon dispersion by 
its leading long-wavelength limit, i?fc ~ ~ t;(fc)|fc|; 




see Eqs. (|2.52| ) and ( |2.53| ). Substituting Eq. ( |3.12| into 
Eq. ( |3.10 ) we finally obtain 



/ he 



cx I \Mcl 



„2 fXl3 



{k,k) ■ ek\\ ■ 
(3.14) 



Formally, this is the leading-order contribution in the 
small parameters v{k)/c\ and hc/{Mc\). Note that the 
magnetic-field dependence of the vertex (3.12) is hidden 



in the canting angle 0, which we approximate by its clas- 
sical value ( 2.36 ). Adding the classical and the hybridiza- 



tion contributions, we finally obtain for the total velocity 
shift, 



(3.15) 



AcA ^ (AeA)o + (Aca)i 
c\ Cx 
Depending on the polarization vector of the phonon 

and on the structure of the first-derivative vector J^^"* 

(2) 

and the second derivative tensor J), , either the first 
or the second contribution on the right-hand side of 
Eq. (3.15) can dominate. In Sec. |IV| we shall explic- 



itly evaluate Eq. (3.15) for the specific parameters of 



(a) 





(b) 



FIG. 5. These Feynman diagrams represent contributions to 
the phonon self-energy which determine the ultrasonic attenu- 
ation rate to leading order in the large-S expansion, (a) One- 
phonon two-magnon processes, (b) Two-phonon one-magnon 
processes. Wavy lines represent the phonon Green's functions 
G'^^°{KX), solid arrows represent the magnon propagators 
Gu-ia^{K), black triangles represent the shifted one-phonon 



two-magnon vertices defined in Eqs. (2.109a 2.109c I, while 
black circles represent the two-phonon one-magnon vertices 
defined in Eqs. (A. 6 A. 7 1. Both processes lead to an at- 
tenuation rate proportional to fe'' at long wavelengths. The 
processes (a) exhibit a stronger dependence on the external 
magnetic field and dominate if the magnon velocities are small 
compared with the phonon velocities. 



of the phonons due to the coupling to the magnon sys- 
tem. In principle, the attenuation rate can be obtained 
using Fermi's golden rule. Due to the rather complicated 
matrix elements of the relevant interaction vertices, we 
find it more convenient to use a many-body approach 
where the damping rate is obtained from the imaginary 
part of the phonon self-energy. The lowest-order inter- 
action processes leading to phonon damping are shown 
diagrammatically in Fig. [5j It turns out that only the 
processes in Fig. [5]ja) involving two magnons in the in- 
termediate lead to a strong magnetic-field dependence of 
the ultrasonic attenuation rate, so that we shall focus 
here on these processes. In the appendix we shall show 
that for CS2CUCI4 where the magnon velocities are small 
compared with the phonon velocities the other processes 
in Fig. [sjb) involving one magnon and one phonon in 
the intermediate state can indeed be neglected. To ob- 
tain the contribution of the processes shown in Fig. [5j^a) 
to the ultrasonic attenuation rate, it is sufficient to re- 
tain only the effective one-phonon two-magnon interac- 
tion process in our coupled magnon-phonon action. We 
thus approximate 



CS2CUCI4 and compare our results with experiments. 



B. Ultrasonic attenuation 

To obtain the magnetic-field dependence of the ultra- 
sonic attenuation rate, we should calculate the damping 



S[X, p, /?] = ^^phoj^j ^ 52mag[/9, P] + ~Sl^^:^[X, /3], 

_ (3.16) 

where the Gaus sian parts S'2p^°[X] and 5'2mag[/3j /?] are 
given in Eqs. (3.3 3.4), while the effective magnon- 
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phonon interaction is 



The diagrams shown in Fig.[5](a) give rise to the following 
contribution to the phonon self-energy, 



1 /A/3t,3t 



k' ■ ^K-K' PkPk' 



2! 



(3.17) 



-ipho 



{K\) 



T 

'n 



E 



If -e^.P 

\^k\k'~k 
2M 

l^k'Jc-k' ' fifcAp 
2M 

+ [K^-K). 



G,„ag(if')Gmag(i^-i^') 



(3.20) 



which corresponds to the effective one- phonon two- 
magnon Hamihonian i?2mag defined in Eq. Kim. Note 
that the vertices in this expression take ttie niagnon- 
phonon hybridization into account. The phonon self- 
energy T,P^^°{KX) and the magnon self-energy Y,tna.g{K) 
are defined in terms of the corresponding Green's func- 
tions: 



Gm.siK) = --{PKpK) = - 



1 



iuj-Ek- i;„iag(i^) ' 



(3.18) 



(3.19) 



Neglecting self-energy corrections to the magnon prop- 
agators, the frequency summation can be easily carried 
out, 



-^pho 



{KX) 



-E 



'^k'%-k ■ ^kxV n{Ek') - n{Ek'-k) 



2M 



iuj — El, 



E, 



k'-k 



I k',k — k' ■ ekxl" niEk') + n{Ek-k') + I 



2M 



-iuj + Ek' + Ek-k 



(3.21) 



where n{E) — l/(e^/"^ — 1) is the Bose function. To 
calculate the phonon damping, we analytically continue 
this expression to the real frequency axis, iw — > a; + iO, 
and take the imaginary part. We obtain 



l^Yt°{u + ^0, fc, A) = -[1 - e-l^]^Y. 



|-p/3+^ „ I 

I-"- k',k'-k ' ^fe^l 

2M 



-6{uj - Ek' + Ek'-k)[l + n{Ek')]n{Ek'-k) 



2M 



- Ek. - Ek-k')[l + n{Ek')][l + n{Ek-k')] 



— [lo ^ — cj, k — fc). 



(3.22) 



The attenuation rate ^kx of a phonon with energy ujkx 
can then be obtained from the imaginary part of the self- 
energy on resonance: 



ImEf°(cjfcA +iO,fc,A) 



IkX 



2lj, 



kX 



At zero temperature we obtain from Eq. (3.22) 



IkX 



1 |f/3*/3t |2 

TT 1 |i k'.k-k' ■ ^kx\ 



2oJkx N 



E 



2M 



X SC^kx — Ek' — Ek-k')- 



(3.23) 



(3.24) 



The leading behavior of ^jkx for small k can be obtained 
analytically. In this limit it is sufficient to use the linear 
approximation (2.52 1 for the magnon dispersion. More- 
over, the interaction vertex f ^ ^, can be approximated 
by the leading long-wavelength limit given in Eq. (2.1101. 



The resulting integration can then be carried out analyt- 
ically and we obtain 



IkX 



64 \2MJ yVBzv^VyJ yi 



Ix{k) 



(3.25) 



' feA 



where rkx = 'v{k)/cx, and the dimensionless function 
/A(fe) is given by 
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2 Sg 



v{k) 



v{k) (2s 



1 - r 



feA 



f^^{k,k)-ek\ +(l-r^A) /a^" (fc, fc±) • e^A 



' feA 



' feA 



fl'^{k) ■ CfcA 

/f^(fc)-efeA f^^{k,k)-e 



kX 



(3.26) 



The vector k± in the last hne is given by kj_ = 
— (vy/vx)kyex + {vx/vy)kxey. We conclude that in the 
regime of small wavevectors considered here, the ultra- 
sonic attenuation rate is proportional to fc"*, which shows 
that in the presence of magnon-phonon interactions the 
phonons remain well-defined quasiparticles. Let us point 
out that in order to obtain the correct fc^-dependence 
of the ultrasonic attenuation rate, it is crucial to take 
the rcnormalization of the one-phonon two-magnon ver- 
tices due to magnon-phonon hybridization into account; 
see Eqs. (2.101a 2.102c). Otherwise one would incor- 



rectly find from Eq. (3.20) that 7fcA reduces to a nonzero 
constant for fc — !■ 0, implying that, due to the coupling 
to the magnon system, the phonons would cease to be 
well-defined quasiparticles. In the appendix, we show 
that the interaction processes represented by Fig. [5jb), 
which involve intermediate states with one magnon and 
one phonon, also lead to the contribution of the order of 
fc^ to the ultrasonic attenuation rate. Fortunately, this 
contribution is negligible in the experimentally relevant 
regime where the magnon velocities are small compared 
with the phonon velocities. 



IV. COMPARISON WITH EXPERIMENTS 

We have measured the longitudinal C22- and C33- 
modes propagating in the plane of the two-dimensional 
layered antiferromagnet CS2CUCI4 down to tempera- 
tures of about 50 mK. Experiments as a function of 
magnetic field were performed in the cone state for 
fields B parallel to the crystallographic a-axis. For 
the experiments presented here, we used a setup which 
allows us to perform simultaneously measurements of 
changes in the ultrasonic velocity Ac/cq and the relative 
attenuation Aa as a function of an external parameter 
like temperature or magnetic field. We employ a 
pulse-echo method using a phase-sensitive detection 
technique.!^ Our measurements can be performed in 
the frequency range of 5-500 MHz. The duration of 
the ultrasonic echo pulse is between 0.1 and 1 /is, while 
the repetition rate is chosen such that it matches the 
available cooling power in the cryostat. It lies between 
100 Hz in the sub-kelvin temperature range, and a few 
kiloherz at higher temperatures. A high-quality single 
crystal of the compound CS2CUCI4 with the size of 



3.68 X 3.8 X 3.65 mm'^ was grown from aqueous solutions 
by an evaporation technique.'^ Two opposite surfaces 
normal to the crystallographic 6-axis and c-axis were 
polished and a pair of piezoelectric thin-film trans- 
ducers was glued to these surfaces. These geometries 
correspond to the longitudinal C22 and C33 acous- 
tic modes, with the wavevector k and polarization e^^ 
parallel to the crystallographic b- and c-axis, respectively. 



To compare the experimental results with our theoret- 
ical predictions, we need realistic estimates for the first 
and second derivatives dJ/dx, dJ' jdr, d'^J/dx^, and 
d'^J'/dr'^ of the exchange couphngs for CS2CUCI4. Recall 
that these derivatives ap pear in the explicit expressions 



(2) 

sor , which in turn enter our theoretical results for 
the elastic constants and the ultrasonic attenuation rate. 
In principle, the spatial dependence of the exchange cou- 
plings can be determined using ah initio methods,'^^'^ 
but quantitative ah initio results for the derivatives of 
the exchange couplings of CS2CUCI4 are currently not 
available. We therefore model the spatial dependence of 
our two relevant exchange couplings J{x) and J'{r) by 
the following phenomenological expressions: 



\2.7l\ and (|2.72a[]2.72c[) for the vector J^^^ and the ten 



J'{r) = j'(d)e-'='('-'"/'' 



(4.1a) 
(4.1b) 



where h and d — /2 are the bond lengths at the 

equilibrium positions, and the dimensionless quantities n 
and k' give the inverse range of the exchange interaction 
in units of the corresponding inverse equilibrium range. 
Since we do not have any a priori knowledge about the 
numerical values of k and k', we determine these param- 
eters by fitting our theoretical predictions for the elastic 
constants to the experimental results. 



According to Eqs. (2.77) and (3.1 ) the classical contri- 



bution (Aca)o to the shift in the phonon velocities aris- 
ing from the motion of the phonons in the classical spin 
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background is determined by the dimensionless ratio 



X 10" 



lim 

\k\^o 



Mcl 



[Jk^ + COS^ Lpo] 



J' 



+Cg [Jk^ cos{Qxb) + —k''^ cos"* lpq cos{Qxh/2)\ 

(4.2) 

for longitudinal phonons in the x-direction (c22-niode), 
and by 



hm ^ 
IfcKo 



pho 



(fc,A) 

2 



-^J'«-sinVo 



X [sl + cos(Q^6/2)] , 



(4.3) 



for longitudinal phonons in the y-direction (cas-mode), 
where J = J(6) and J' = J'{d). Moreover, to evalu- 
ate the hybridization contribution (Aca)i to the velocity 
shift given in Eq. (3.14) and the ultrasonic attenuation 
rate given in Eqs. (3.25 3.26) we need the dimension- 
less scalar products /^(fc) • ek\ and /^^(fc,fc) • CfcA, 
where the vector functions f^^ {k ) and f2^{k,k') are 



defined in Eqs. (3.13a) and (3.13b|. To explain the data 



for the low-temperature measurements of the longitudi- 
nal C22 and C33 phonon modes, we set k = [k^, ky) — e^x- 
In this case the relevant scalar products simplify to 



kl [Jn{l - cos{Qxb)) 

+ 2J'k' cos^ .^0(1 - cos(Q^6/2))] 
+ kpj'ti' sin^ (^0(1 - cos(Q^6/2))| , 

f^\k,k)-k^^[klh[jKsMQ.,h) 
v[k) ^ 

+ J'k' snY{Qxb/2) cos^ (p\ 
-f 3fc2;fc^cJ'K' sin^jQ cos (/3o sin(Qa;5/2)| , (4.5) 



(4.4) 



f2^{k,k^_) 



w(fc) I 



— kyh[jKSiT\{Qxh) + J'k! cos^ (^q sin((5^6/2)] 

X cJ'k' sin (/7o cos iy9o sin((5 a; &/2) I , (4.6) 

where we have used the fact that the wavevector Q of 
the spiral has only an a;-component. 

Setting k = (1,0) for the C22-niode and k — (0,1) 
for t he C33 -mode, we can explicitly evaluate the terms in 
Eq. (3.15) and calculate the velocity shifts. In Fig. [6] we 
compare our theoretical results for the magnetic-field de- 
pendence of the velocity shifts of the longitudinal phonon 
modes with our experimental data at low temperature. 
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FIG. 6. (Color online) Measured velocity shifts of the longitu- 
dinal C22-phonon mode (squares) and the css-mode (circles) of 
CS2CUCI4 as a function of the magnetic induction B taken at 
T = 52 mK (c22-mode) and T = 48 mK (cas-mode). The solid 
line is a fit of the data for the C22-mode to our theoretical re- 
sults at zero temperature given in Eqs. (2.77 |3.1 3.14 



3.151 



From this fit we obtain for our phen omen ologic al parame- 
ters K and k! introduced in Eqs. (4. la I and (4. lb I the values 
K ~ 15 and hi ~ 51. Using the same values for the C33-mode, 
we obtain the theoretical prediction represented by the dashed 
line. From the inset it is obvious that in the weak-field regime 
B < 2.5 Tesla where our spin-wave approach should be most 
accurate we obtain excellent agreement between theory and 
experiment. The anomalies near the critical field Be ~ 8.5 
Tesla might have a different physical origin than magnon ex- 
citations, so that these features cannot be explained within 
our spin-wave approach. 



The data for the C22-inode was obtained at T = 52 mK, 
while the data for the C33-mode was taken at T = 48 mK. 
The parameters k and k! were determined by fitting our 
theoretical results for the C22-nrode to the data. The re- 
sulting values K w 15 and k' w 51 are then inserted back 
into our expression for the C33-mode, so that our theo- 
retical prediction for the C33-mode does not contain any 
adjustable parameters.'^ From the inset in Fig. |6] it is 
obvious that in the weak- field limit B <2.b Tesla, where 
our calculations of the velocity shifts are expected to be 
most accurate, we obtain excellent agreement between 
theory and experiment. The deviations between the ex- 
perimental data and our calculations for larger fields sig- 
nal the breakdown of our theoretical approach which does 
not take into account higher order fluctuation corrections 
and other types of excitations. These are likely to play 
a role in the vicinity of the critical magnetic field where 
the magnetic order vanishes. Near the critical field for 
temperatures below 0.3 K the ultrasonic attenuation ex- 
hibits a double peak structure which will be discussed in 
a separate publication. 

Having fixed the fit parameters k and k' from the ve- 
locity shifts, our theoretical result for the ultrasonic at- 



tenuation rate given in Eqs. (3.25) and (3.261 does not 
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FIG. 7. (Color online) Experimental results for the relative 
ultrasonic attenuation Aa in CS2CUCI4 of the longitudinal 
C22-phonon mode (squares) and the csa-mode (circles) taken 
at T = 52 mK (c22-mode) and T = 48 mK (css-mode). The 
solid and the dashed lines represent our corresponding theo- 
retical predictions given in Eqs. (3.251 and ( |3.26[ |. The atten- 
uation rates have been obtained using the same values k = 15 
and k' = 51 for the fit parameters as in Fig. [g] The constant 
offset in the attenuation rates is due to the fact that with our 
experimental technique we can only measure relative attenua- 
tion rates. Note that for very small fields the data approaches 
Aa = 0. 



contain any adjustable parameters. In Fig. [T] we com- 
pare the results of our calculations for the attenuation 
rate of the C22 and C33 phonon modes with the experi- 
mental data of the relative attenuation Aa. For fields 
in the range h < O.Shc the measured attenuation rates 
are rather small, while one observes a strong enhance- 
ment for h ^ he- This indicates that, at least for small 
magnetic fields, other effects leading to the large atten- 
uation background are not field dependent. The overall 
shape of the data is reproduced by our theoretical curves 
rather well, although our spin-wave calculation cannot 
reproduce the shoulder just below the critical field. Note 
that according to Eqs. (3.25) and (3.26) the attenuation 



rate for fields in the vicinity of the critical field can be 
approximated by 



7fcA 



64 \2m) 



(4.7) 



where we have approximated sin « 1 and cos 6 w 
yj2{l — h/hc) for h close to he, and have used the fact 
that in CS2CUCI4 the ratio r^x ~ v{k)/c\ is small com- 
pared with unity. One should keep in mind, however, 
that in the vicinity of higher order fluctuation cor- 
rections which we have neglected within our spin-wave 
expansions are likely to become important, so that we do 
not expect that Eq. ( |4.7[ ) is quantitatively accurate very 
close to he- Nevertheless, from Fig. [7| it is clear that the 
factor of (1 — h/hc)~'^ in Eq. (4.7) gives a satisfactory 



description of the strong enhancement of the attenuation 
rate in the vicinity of the critical field. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have presented both theoretical and 
experimental results for the magnetic-field dependence of 
the elastic constants and the ultrasonic attenuation rate 
in the cone state of the frustrated quantum antiferro- 
magnet CS2CUCI4. Our calculations are based on the as- 
sumption that two-dimensional magnons are well defined 
excitations in this material, so that we may use the 1/S 
spin- wave expansion in combination with an expansion in 
powers of phonon operators. Using a simple phenomeno- 
logical parametrization of the spatial variations of the 
exchange couplings involving two adjustable parameters, 
our theoretical results for the magnetic-field dependence 
of the elastic constants agree quite well with our ex- 
perimental measurements, in particular in the weak-field 
regime B < 2.5 Tesla where our perturbative approach is 
expected to be most accurate. Our theoretical results for 
the ultrasonic attenuation rate reproduce the strong en- 
hancement close to the critical magnetic field, although 
our approach is likely to break down in the vicinity of 
the quantum critical point. 

The fact that the ultrasonic attenuation rate is pro- 
portional to implies that the coupling to the magnon 
excitations does not destroy the quasiparticle charac- 
ter of the long- wavelength acoustic phonons. To arrive 
at this result, it was essential to work out the relevant 
magnon-phonon matrix elements with great care, taking 
into account that due to magnon-phonon hybridization 
in the cone state the cubic part the pure magnon Hamil- 
tonian indirectly contributes to the effective one-phonon 
two-magnon interaction. Up to now the predicted fc**- 
dependence of the attenuation rate on the phonon mo- 
mentum could not be detected experimentally. A possi- 
ble reason is that the magnetic contribution is small com- 
pared with the background, which does not depend on 
the magnetic field. In the vicinity of the critical field, the 
magnetic contribution becomes large, but it is question- 
able whether our calculation is still valid in this regime. 

We have intentionally restricted the calculations pre- 
sented in this paper to the magnetically ordered cone 
state of CS2CUCI4, because in this case our spin- wave 
approach is well justified. The investigation of magnon- 
phonon interactions in the immediate vicinity of the 
quantum critical point and in the other phases of 
CS2CUCI4 such as the spin-liquid phase or the ferromag- 
netic phase is left for future work. Since magnons in 
frustrated magnet s can have anomalous properties such 
as strong damping,'2SllSI to magnon-magnon interac- 
tions, it should be interesting to take a closer look at 
these interaction processes between magnons in the cone 
state, especially in the vicinity of the quantum critical 
point. 
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(A.6) 
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To second order in the two-phonon one-magnon vertex, 
we obtain for the phonon self-energy 



APPENDIX: PROCESSES INVOLVING ONE 
PHONON AND ONE MAGNON IN THE 
INTERMEDIATE STATE 

In this appendix, we show that in CS2CUCI4 where the 
phonon velocities are large compared with the magnon 
velocities the ultrasonic attenuation rate due to processes 
involving one phonon and one magnon in the intermedi- 
ate state shown in Fig. [sjb) is small compared with the 
contribution from processes with two magnons in the in- 
termediate state shown in Fig.jsjja). 

Consider the part iJ^^^g the magnon-phonon inter- 
action involving two powers of the phonon operators and 
one power of the magnon oper ators. It can be obtained 
from given in Eq. (2.601 by setting n — 2 and ex- 

panding the spin operators to first order in the Holstein- 
Primakoff magnon operators. The Fourier transform of 



(2) 

the relevant phonon part C/^ ^, is given in Eq. (2.66). We 
obtain 
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2pho 
Imag 



^k\,k'\'^k\Xk'X'b^k+k 
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where ^^xk'X' ~ ^kx^kk/'^^'^'X'- The Matsubara sums 
can now be carried, out and we obtain 



^pho 



t 1^ kX-k'X' I 

N^, 2ujk'X'M^ 

k' X' 



n{LL>k' 



n{Ek-k') 
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n{uJk'X') + n{Ek-k') + 1 
iuj — LOk'X' — Ek-k' 



iuj + LUk'x - Ek-k' 
[K^-K), (A.9) 



At zero temperature this yields for the ultrasonic atten- 
uation rate 



vr 1 



IkX 



^xxp 

kX.-k'X' 



2ukxN^, 2uk'X'M^ 



-^{^k\ — ^k'X' — Ek-k')- 

(A.IO) 



In the long-wavelength limit we obtain for the relevant 
matrix element of the interaction vertex 



where the vertices are 

j^XXb^ 
^ kX.k'X' 

T^XXb _ I -pAAD 
^ kX,k'X' — fifcA-"- kk' ^k'X 

P give 

the direction labels. 



^kX^ kk' fife' A' 
,t T^XXb, 



(A.2) 
(A.3) 



with T^^^ and T-^^'^ given by the following matrices in 



k+k',Q 



r(2-) 
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t(2+) _ 
k,Q 



t(2-) 
k,Q 

,(2+) , 



fe'.Q 



l(2^ 



j^XXb 
'-kk' y 



(A.4) 



where the tensors J^. q were introduced in Eqs. (2.78b 
2.78a). From the symmetries in Eq. (2.79) it follows that 



for small k and k' the vertices F^^^ and T^^^ vanish 
quadratically in k and k' . Transforming to the Bogoli- 
ubov basis, our Hamiltonian ( |A.1[ ) assumes the form 

^Imag = 77-7^ E E [^kX^k'X'^kxXk'X'l3i+k' 

rfcA^fe^A'^-feAA:_fc'A'/3fc+fc' , (A. 5) 
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[(fc-fe')-VQ](fe-VQ)(fc'-VQ)J^^ 



(A.ll) 



To carry out the integration in Eq. ( A.IO ), we use the fact 



that in the experimentally relevant regime the phonon ve- 
locities are much larger than the magnon velocities. Then 
the scattering surface defined by uJkx ~^k'X' — Ek-k' = 
can be explicitly calculated to leading order in v{k)/c\. 
Using circular coordinates we obtain the parametric rep- 
resentation 



cx' 



(A.12) 



where k'{ip') — cosLp'sx + sin(^'e.y, and 



u{k, ip') = \l v1{kx — cos ^'Y -\- v^{ky — sin^j')^. 



(A.13) 
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To leading order in v{k)/c\ we obtain 



7fcA 



4 y2Mj \VbzJ ^ \Mci,J VcA 



'^V^^^l4Ar^^^(fc,^')e.'Vp, (A.14) 



where the dimensionless matrix element is given by 



'Q 



2u{k,ip') ca' ^ 



(2) 



(A.15) 



r 



Comparing Eq. (A.14) with the corresponding damping velocities, the processes with two magnons in the inter- 
rate (3.25) due to one-phonon two-magnon processes, we mediate states are dominant, although both processes 
see that in the experimentally relevant regime where the yield contributions to the ultrasonic attenuation rate 
magnon velocities are small compared with the phonon which are proportional to k^ at long wavelengths. 



^ For a collection of recent reviews see, for example, 
C. Lacroix, P. Mendels, and F. Mila, Introduction to Frus- 
trated Magnetism, Springer Series in Solid-State Sciences 
164 (Springer, Berlin, 2011). 
^ A. N. Nasyrov, H. Shodiev, Z. Tylczynski, A. D. Karaev, 

and V. S. Kim, Ferroelectrics 158, 93 (1994) 
^ R. Coldea, D. A. Tennant, R. A. Cowley, D. F. McMorrow, 
B. Dorner, and Z. Tylczynski, J. Phys.: Condens. Matter] 
[~ 8, 7473 (1996)j 
" R. Coldea, D. A. Tennant, R. A. Cowley, D. F. McMorrow, 
B. Dorner, and Z. Tylczynski, |Phys. Rev. Lett. 79, TsT] 



(1997) 

S. Bailleul, J. Holsa, and P. Porcher, Eur. J. Solid State 
Inorg. Chem. 31, 431 (1994). 

R. Coldea, D. A. Tennant, K. Habicht, P. Smeibidl, 
C. Wolters, and Z. Tylczynski, Phys. Rev. Lett. 88, 137203| 
(2002) 

T. Radu, H. Wilhelm, V. Yushankhai, D. Kovrizhin, 
R. Coldea, Z. Tylczynski, T. Liihmann, and F. Steglich, 



[Phys. Rev. Lett. 95, 127202 (2005)t 
^ M. Y. Veillette, J. T. Chalker, and R. Coldea, |Phys. Rev. 



B 71, 214426 (2005) 
" M. Y. Veillette, A. J. A. James, 



and F. H. L. Essler, Phys. 



Rev. B 72, 134429 (2005) 

^" M. Y. V eillette and J. T. Chalker, |Phys. Rev. B 74 , 052402 



"(2006)1 



D. Dalidovich, R. Sknepnek, A. J. Berlinsky, J. Zhang 
and C. Kallin, 'Phys. Rev. B 73, 184403 (2006) 



D. L. Kovrizhin, V. Yushankhai, and L. Siurakshina, |Phys. | 
Rev. B 74, 134417 (2006) 

M. Q. Weng, D. N. Sheng, Z. Y. Weng, and R. J. BursiU, 



Phys. Rev. B 74, 012407 (2006) , 
^'^ M. Kohno, O. A. Starykh, and L. Balents, [Nat. Phys. 3, 



790 (2007) 

O. A. Starykh and L. Balents, |Phys. Rev. Lett. 98, 077205 



(2007), _ 

M. Kohno, Phys. Rev. Lett. 103, 197203 (2009) 



O. A. Starykh, H. Katsura, and L. Balents, Phys. Rev. B 
82, 014421 (2010)1 



C. H. 



Phys. 



Chung, J. B. Marston, and R. H. McKenzie, |X] 

Condens. Matte r 13, 5159 (2001)., 

J. O. Fjaerestad, W. Zh eng, R. R. P. Singh, R. H. Mc Ken- 
zie, and R. Coldea, Phys. Rev. B 75, 174447 (2007) , 
B. Liithi, Physical Acoustics in the Solid State (Springer, 
Berhn, 2005). 

A. Sytcheva, O. Chiatti, J. Wosnitza, S. Zherlitsyn, A. A. 
Zvya gin, R. Coldea, and Z. Tylczynski, |Phys. Rev. B 80, | 



224414 (2009) 

Tachiki and S 



M. 



(1974) 



Maekawa, Prog. Theor. Phys. 51, 
K. Jung, 



SoUd State Commun. 132, 393 



^ATR. Lim and J 
(2004) 

S. Aplesnin, IPhys. Lett. A 333, 446 (2004) 
J. H. Kim and J. H. Han, Phys. Rev. B 76, 054431 (2007 
^® Y. Zhou and P. A. Lee, |Phys . Rev. Lett. 106, 056402 



(2011) 



R. Coldea, D. A. Tennant, and Z. Tylczynski, Phys. Rev. 
B 68, 134424 (20 03) 

F. Schiitz, M. Ko Uar, and P. Kopietz, |Phys. Rev. Lett. 



I 91, 017205 (2003) 

^I. Spremo, F. Schiitz, P. Kopietz, V. Pashchenko, B. Wolf, 
M. Lang, J. W. Bats, C. Hu, and M. U. Schmidt, . Phys.| 
[~ Rev. B 72, 174429 (2005) 

K. Foyevtsova (private communication) 
N. Hasselmann and P. Kopietz, Europhys. Lett. 74, 1067 



(200 6) 
A. Kreisel, F. 



Sauli, N. Hasselmann, and P. Kopietz, |Phys. 



Rev. B 78, 035127 (2008) 

A. Perelomov, Generalized Coherent States and Their Ap- 
plications (Springer, Berlin, 1986). 

B. Liithi, G. Bruls, P. Thalmeier, B. Wolf, D. Finsterbusch, 
and I. Kouroudis, ,J. Lo w Temp. Phys. 95 , 257 (1994) 

N. Kriiger, S. Belz, F. Schossau, A. A. Haghighirad, P. T. 
Cong, B. Wolf, S. Gottlieb-Schoenmeyer, F. Ritter, and 

W. Assmus, Cryst. Growth Des. 10, 4456 (2010) 

Y.-Z. Zhang, H. O. Jeschke, ' 
78, 205104 (2008) 
K. Foyevtsova, I. Opahle, Y. Zhang, H. O. Jeschke, and 



and R. Valentf, |Phys. Rev. B 



20 



R. Valentf, Phys. Rev. B 83, 125126 (2011). 
^® Note that the cas-phonons do not modulate the bond 
lengths of the coupHng J, so that the associated velocity 
shift depends only on k . 



^® A. L. Chernyshev and M. E. Zhitomirsky, Phys. Rev. Lett. 
97, 207202 (2006) 



A. L. Chernyshev and M. E. Zhitomirsky, Phys. Rev. B 



79, 144416 (2009)| 



